CIRRUS PARCEL MODEL COMPARISON PROJECT PHASE 1 

Ruei-Fong Lin 1 * David O’C. Starr 2 , Paul J. DeMott 3 
Richard Cotton 4 , Eric Jensen 5 , and Kenneth Sassen 6 

Universities Space Research Association, Greenbelt, MD, USA 
2 Laboratory of Atmosphere, NASA/GSFC, Greenbelt, MD, USA 
3 Dept. of Atmospheric Sciences, Colorado State University, Fort Collins, CO, USA 
4 Meteorological Research Flight, Hants, UK 
5 N AS A/ Ames Research Center, Moffett Field, CA, USA 
6 Dept, of Meteorology, University of Utah, UT, USA 


1 INTRODUCTION 

The Cirrus Parcel Model Comparison (CPMC) is a 
project of the GEWEX Cloud System Study Work- 
ing Group on Cirrus Cloud Systems (GCSS WG2). 
The primary goal of this project is to identify cirrus 
model sensitivities to the state of our knowledge of 
nucleation and microphysics. Furthermore, the com- 
mon ground of the findings may provide guidelines 
for models with simpler cirrus microphysics modules. 


Table 1: Simulation identifiers. 
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We focus on the nucleation regimes of the warm 
(parcel starting at — 40° C and 340 hPa) and cold 
(-60°C and 170 hPa) cases studied in the GCSS 
WG2 Idealized Cirrus Model Comparison Project 
[Starr' et al, 2000]. Nucleation and ice crystal 
growth were forced through an externally imposed 
rate of lift and consequent adiabatic cooling (Ta- 
ble 1). The background haze particles arc as- 
sumed to be lognormally-distributed H 2 SO 4 parti- 
cles. Only the homogeneous nucleation mode is al- 
lowed to form ice crystals in the HN-ONLY runs; 
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all nucleation modes are switched on in the ALL- 
MODE runs. Participants were asked to run the 
HN-A-fixed runs by setting A = 2 (A is further dis- 
cussed in section 2) or tailoring the nucleation rate 
calculation in agreement with A = 2 1 . The depth 
of parcel lift (800 m) was set to assure that parcels 
underwent complete transition through the nucle- 
ation regime to a stage of approximate equilibrium 
between ice mass growth and vapor supplied by the 
specified updrafts. 

2 MODEL DESCRIPTIONS 

Five parcel modeling groups participated in the 
CPMC (Table 2). Hereafter, we will refer to these 
models as the C, D, J, L, and S models, respectively, 
as denoted in the table. 

The estimate of the nucleation rate of ice in solu- 
tion droplets, Jhaze 1 remains an active research area. 
Jhaze was computed using either (1) the modified 
classical theory approach (model J) or (2) the effec- 
tive freezing temperature approach (hereafter, T e fj 
models, models C, D, L, S). 

The T e fj models attempt to directly link mea- 
sured Jhaze to nucleation rates of equivalent-sized 
pure water droplets J w via the effective freezing tem- 
perature, which is defined as 

r eff = T + AAT m , (1) 

such that Jhaze = Jw(T e ff) as introduced by Sassen 
and Dodd [1988]. In (1), A T m is the equilib- 
rium melting point depression (positive valued), 
which depends on solute wt%, and A is an empir- 
ical coefficient to account for additional suppres- 
sion/enhancement of nucleation temperature due to 

l Note that A = 2 agrees approximately with data pre- 
sented by Koop et al. [1998]. 
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Table 2: Participant cirrus parcel models. 


Organization 

Investigator 

UKMO 
Cotton (C) 

CSV 

DeMott (D) 

ARC 

Jensen (J) 

GSFC 
Lin (L) 

U. Utah 
Sassen (S) 


Bin characteristic" 

discrete 

continuous 

continuous 

continuous 

particle tracing 


Haze size 6 

r e<? or % 

T eq 

V eq 

1 eq 

or lit 


A 

1.5 

1.5 

varying" 

1.0 

1.7 


deposition coef. /3* 

0.24 

0.04 

1 

0.1 

0.36 


References 

Spice et al. 
[1999] 

DeMott 
et al. [1994] 
DeMott 
et al. [1998] 

Jensen et al 
[1994] 
Tabazadeh 
et al [1998] 

Lin [1997] 

Sassen and 
Dodd [1988] 
Sassen and 
Benson [2000] 

.ill * 



size bin or a certain distribution of particle sizes is allowed in a bin. 

6 r e , vs. £ denotes either using the equlibrium-sized haze approximation or computing the diffusional 
growth of haze particles explicitly. 
c See section 2 for detailed discussion. 


non-ideal interaction between ions and condensed 
water. Although Sassen and Dodd [1988] noted that 
an average A for different solutions was around 1.7, 
values for specific solutions may range from 1 to 2.5. 



Temperature (C) 


Figure 1: JhazeV vs. temperature for solute wt,% 
5, 15 and 25%. Solid, dashed, dash-dotted, dotted 
curves denote models J, C, S, models D and L (same 
curves), respectively, for A — 2. 

In model J, recent direct data on ice/solution sur- 
face tension was incorporated and activation energy 
was inferred from recent laboratory measurements 
of Jhaze for H 2 SO 4 particles following Tabazadeh et 
al [1997] and Koop et al [1998]. This approach to 
determine Jhaze can be interpreted as a T e // scheme 
with varying A (Figure 1). The intrinsic A varies in- 
versely with solute wt% and temperature. Also, the 
differences in the sensitivity of JhazeV (V is the vol- 
ume of the particle) to solute wt% between these two 
approaches may lead to systematic differences in the 


freezing haze size distributions. Nucleation rate data 
over a wide range of values, e.g., data points beyond 
critical freezing conditions, are needed to diminish 
the inconsistency between the two approaches. 

Little constraint was imposed on formulating het- 
erogeneous nucleation because theoretical and ex- 
perimental understanding are still quite poor. Mod- 
els C and L employ ice saturation ratio dependent 
parameterizations of activated IN following Spice et 
al [1999] and Meyers et al [1992], respectively. 
These parameterizations are expected to represent a 
maximum heterogeneous nucleation impact. 

Haze particles of the given H 2 SO 4 aerosol distri- 
bution are subject simultaneously to heterogeneous 
and homogeneous nucleation in models D and S. The 
number concentration of the activated IN in model D 
is computed following DeMott et al [1998] based on 
field experiment data. This treatment was expected 
to yield the most conservative estimate of IN in cir- 
rus. Model S computes the activated freezing nuclei 
using T e ff dependent Fletcher equation [5assen and 
Benson , 2000], where parameters were set to yield 
the most favorable conditions for heterogeneous nu- 
cleation. 

Participants either assumed that haze particles 
are in equilibrium with the environment or com- 
puted the diffusional growth of haze particles di- 
rectly (Table 2). The diffusional growth rate of haze 
particles more or less exponentially decreases with 
temperature as caused by water vapor saturation 
pressure. The response time scale to the deviation 
from equilibrium can be considerably greater than 
one model time step in a swift updraft in a cold 
environment. Therefore, large haze particles may 




become more concentrated than the corresponding 
equilibrium-size particles in such conditions. This 
may result in considerable delaying of haze growth 
in models C and S (Table 2) and affect ice particle 
formation rate. 

3 RESULTS AND DISCUSSIONS 

As we proceed to describe the results and differences 
between models, it must be noted that the bench- 
mark is not necessarily the median or the average of 
model results. The predicted N t (ice number concen- 
tration) at 800 m above the starting point is com- 
pared (Fig. 2). In the HN-ONLY cases, to a first 
order approximation, the logarithm of N{ increases 
quasi-linearly with the logarithm of updraft speed. 
The predicted N t by models D, S and L are close; 
N x by models J and C form the lowest and highest 
bounds in the six cases, respectively. 
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Figure 2: N{ predicted vs imposed updraft speed. 
The unfilled and filled bars denote HN-ONLY and 
ALL-MODE, respectively. 


Cirrus initiation occurred over a narrower range 
of altitude and RH X (relative humidity over ice) in 
the warm HN-ONLY cases than in the cold cases 
(Fig. 3). The increasing sensitivity of the cloud base 
RH t as temperature decreases in the four T e ff mod- 
els is primarily caused by A. 

Heterogeneous nucleation is a possible explana- 
tion of the discrepancy between the observed thresh- 
old RH W for cirrus formation and the theoreti- 
cally derived threshold RH W (relative humidity over 
water) for homogeneous nucleation of H 2 SO 4 or 
(N H4) 2SO4 solution particles; e.g., [Heymsfield and 
Miloshevich , 1995]. Cirrus properties are affected 
by the dominant nucleation mode in cloud initia- 
tion because of the distinct characteristics of the two 
modes. 

The cloud base height, RH X and peak RH t in the 
ALL- MODE cases (not shown) vary even more be- 
cause of our respective unbounded choices of hetero- 
geneous nucleation. The impact of heterogeneous 
nucleation on lowering N t , peak RH X , and cloud for- 
mation altitude is extremely sensitive to the onset 
conditions for nucleation and the subsequent ice par- 
ticle formation rate. With heterogeneous nucleation, 
the peak RH t is lower in all but the case WalOO by 
model S. The predicted N t is reduced in all but the 
case Ca004 by model S. 

We now discuss the results of the HN-A-fixed sim- 
ulations. The nucleation regimes of Wh020L and 
Ch020L take place within the temperature range of 
-43.2 to -44.2°C and -63.2 to -64.2°C, respectively. 
The effect of temperature variation on nucleation 
rates within this 1°C range is secondary, compared 
to the evolution of haze solute wt%. Thus, it is jus- 
tified to analyze and visualize results according to 
the z — z\y coordinate (Fig. 4). 



Figure 3: The RH t at cloud base Z& ( Ni = 1 liter -1 ) 
and the corresponding A RH t , defined as the differ- 
ence between peak RHi and RH t at zt (the HN- 
ONLY cases). 


Figure 4: Ice water content (IWC), N x , ice particle 
formation rate and RH W as functions of z - z\>. 

The triggering RH W range was reduced signifi- 
cantly, to less than 2% in Wh020L and 5% in Ch020L 





in comparison to 3% and 8% in Wh020 and Ch020. 
The predicted TV* is only marginally affected. 

At the beginning of the nucleation stage in 
Wh020L, ice particle formation rates by the four 
T e fj models are close. However, models C and D 
reach much larger RH W that leads to larger instan- 
taneous nucleation rates, and maintain the peak ice 
formation rate longer than the other two models. 

Quite contrarily, the Ni curves of models D and L 
in Ch020L distinctly separate from those of models 
C and S. This grouping incidentally coincides with 
the grouping according to the haze size specifica- 
tions. Large haze particles are more concentrated 
than the corresponding equilibrium values in mod- 
els C and S. Yet, the nucleation regime in model S 
was not sustained as long as in model C; a similar 
finding is noted when comparing results of model D 
and L. The results of model J feature slow ice par- 
ticle formation rate, long nucleation duration, and 
broader freezing haze number distribution. 

The above results indicate that nucleation dura- 
tion time and the maximum nucleation rate achieved 
are the two key components in determining the final 
Ni. These two factors are sensitive to the growth 
rates of small ice crystals, which under the influence 
of the kinetic effect are sensitive to the deposition 
coefficient, f3 t . It was found that varying (i x from 
0.04 to 1 (Table 2) would result in about a factor of 
4^5 (Wh020L) and 9~12 (Ch020L) variation in N % 
by models C and L. 

4 SUMMARY 

Results of Phase 1 of CPMC projects show that the 
predicted cloud properties strongly depend on up- 
draft speed. Significant differences are found in the 
predicted N x . Detailed examination revealed that 
the homogeneous nucleation formulation, aerosol 
size specification, ice crystal growth (especially the 
specification of the deposition coefficient for ice) and 
water vapor uptake rate were the critical compo- 
nents. These results highlight the need for new lab- 
oratory and field measurements to infer the correct 
values for critical quantities in the cirrus regime. 

No attempt was made to scrutinize the causes of 
differences in ALL-MODE simulations due to the 
substantial differences in formulation of heteroge- 
neous nucleation. Nevertheless, it was confirmed 
that the expected effect of a heterogeneous nucle- 
ation process is to decrease N % and the RH l required 
for cloud initiation. Clearly, new measurements of 
ice nuclei activation in cirrus conditions are war- 
ranted. 

CPMC Phase 1 was conducted based on a single 


CCN distribution. Phase 2 of the CPMC, now un- 
derway, examines the effects of varying aerosol dis- 
tributions. Sensitivity of model results to CCN com- 
position is indirectly made by altering A. 
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